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Abstract. The generalized second-order partial derivatives of where r is 
the radial distance in 3D, are obtained using a result of the potential theory of 
classical analysis. Some non-spherical-regularization alternatives to the standard 
spherical-rcgularization expression for the derivatives are derived. The utility of 
a spheroidal-regularization expression is illustrated on an example from classical 
electrodynamics. 



1. Introduction 

The expression for the Laplacian of 1/r, 

V 2 - = -4ac6(r\ (1) 
r 

where r — \r\ = (x 2 + x\ + x 2 ) 1 ^ 2 is the magnitude of a vector r = (xx, X2, £3), 
the Laplacian V 2 is the differential operator d 2 jdx\ + d 2 /dx 2 + d 2 /dx 2 and S(r) = 
S(xi)8(x2)S(x3) is the 3D delta function, is well-known and its correct use involves 
only the elementary rules of the delta-function formalism. The Laplacian (JIJ is useful 
in innumerable calculations of electrodynamics, but sometimes the 2nd-order partial 
derivatives (d 2 / dxidxj){l/r) themselves are needed. Examples are the calculation of 
the electromagnetic fields of point dipoles [T] and solving the Poisson equation for the 
difference between the Coulomb- and Lorenz-gauge vector potentials of a uniformly 
moving point charge [2]- 

While the expression for (d 2 /dxidxj)(l/r) has been known for some time, its 
correct use is rather more intricate than that of V 2 (l/r). This expression is usually 
written as [3] 

d 2 1 3x i x 1 ■- r 2 5 l1 4?r 

— 8(r)8ij, (2) 



dxidxj r r 5 3 



which hides an important fact that an integration of the product of the first term on 
the right-hand side and a 'well-behaved' test function over a 3-dimensional domain 
that includes the origin r = still has to be regularized because of the term's 1/r 3 
behaviour at the origirjj The regularization assumed in ([2J is of a specific, 'spherical' 



\ An expression for the (generalized) 2nd-order derivative of 1/r that takes this circumstance into 
account explicitly seems to have been given first in |4], p 28. 
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kind, which can be effected in many equivalent forms, e.g. 

3xjX^ r Sij 3x^Xj v S^j , . 

J — J - ->• whm — — , J (3) 

r 5 e^o (r 2 + e 2 ) 5 / 2 

-» wlim 3XtXj T ^ 6(r-e). (4) 

Here, the symbol wlim indicates the weak limil[§| and O(-) is the Heaviside step function. 
The form Q is implemented automatically when spherical coordinates are used in 
integration and the angular integration is done first. 

Relation ([2]) was derived as a delta- function identity by Frahm [3j. In the 
present paper, we obtain it easily as a generalized (distributional) derivative using 
a result of the potential theory of classical analysis. We also derive some 'non- 
spherical' alternatives to @; while integration in spherical coordinates provides a 
straightforward implementation of the spherical regularization implied in ([5]), the 
use of a non-spherical regularization may be more advantageous computationally in 
some applications. An example from classical electrodynamics where a spheroidal 
regularization is useful is given in Appendix. Apart from the utility, awareness of 
non-spherical alternatives is important for the correct use of the standard, 'spherical' 
expression. 

Our treatment assumes no knowledge of the theory of generalized functions 
and generalized (distributional) derivatives beyond the elementary delta-function 
formalism, but it should help elucidate the operational meaning of some essentially 
non-classical mathematical objectt[jj|. Presenting the 2nd-order partial derivatives of 
1/r from the very beginning as generalized (distributional) derivatives should help 
avoid the pitfalls that would await anyone attempting to use the expression ^ , or its 
very recent 'generalization' 6 (remarked on in the last section), in an integration in 
non-spherical coordinates. This topic and its treatment are suitable for graduate and 
advanced undergraduate courses of electrodynamics. 

A generalized-function treatment of the singularities that may arise at the origin 
r = of a spherical (polar) coordinate system has been presented recently in this 
Journal by Gsponer [7]. In contrast to Gsponcr, our approach does not rely on any 
specific choice of coordinates. 



2. Derivation using derivatives of the potential of an extended density 

Let 

</>(r) = J d3r '^Jf-> R = \r-r'\ (5) 

be the 'potential' created by a density p(r) that is assumed to be a 'well-behaved' 
localized function of r. While the 2nd-order partial derivatives d 2 <f>(r) / dxidxj 
perfectly exist at any point r if the density p(r) is sufficiently 'smoothjjl they cannot 
be calculated by a straightforward differentiating inside the integral that defines the 

§ wUm a ->a 9a{r) = g(r) iff lirn a ->a / f(r)g a (r) = f d 3 r f(r)g(r), where g(r) is in general a 
generalized function (distribution) and f(r) is any well-behaved test function. 

|| For a general-function theoretic underpinning of our approach, the interested reader is directed 
to the canonical regularization method of Gelfand and Shilov and the completeness theorem of 
generalized functions [5]. 

If the density p a (r ) has a sharp surface, the 2nd-order derivatives of the potential are discontinuous 
at that surface. 
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potential because the resulting integrand is not integrable at the point r' = r. A 
correct way of performing here the differentiation under the integral sign is given by 
the following formula: 

d 2 <t>(r) _ Hm f ^ 8 2 1 4tt 

dxidxj ^oJ R>e dxidxj R 3 lJ 

/i , , 3(xi — x'Alx-j — x'A — R 2 5a Air 
dV p(r')^ J R5 3 > ^ Q(R - e) - - p(r)«%. (6) 

Here, the integration domain excludes a ball of radius e and centre at r' = r; the 
2nd line is a more explicit transcription of the 1st line, with the integration domain 
R = \r' — r\ > e expressed by the Heaviside step- function factor Q(R — e) in the 
integrand. An instructive derivation of formula ([6]) is given in the well-known text 
of Tikhonov and Samarskii [8] (they consider the most involved case i = j, but 
the generalization to any i,j = 1,2,3 is straightforward). This derivation uses the 
divergence theorem with a function that has only a 1 / R singularity at the point R = 
and thus, unlike the informal proof in [3J, is fully legitimate in classical analysisQ. 

Note that the regularization of the integral on the right-hand side of © is 
explicitly of the spherical kind. If the excluded integration domain was assumed 
to have a non-spherical shape, the ensuing regularization would in general yield a 
different value for the integral, resulting in an incorrect value for the derivative in 
question. The reason why different regularizations yield in general different values of 
the regularized integral is that, because of the factor (d 2 /dxidxj)(l/R), the integrand 
goes through large positive and negative variations near the point R = 0. The integral 
thus can be made convergent only conditionally [8]. 

Defining a generalized (distributional) derivative 

B 2 1 3(x t ~x' i )(x j -x'A-R 2 6 lJ 4tt 



a**, i s *s - — «**-«>- T '(- -')% m 

(denoted by a bar to distinguish it from a classical derivative [H UHl E] ) > formula © 
can be written simply as 

dxtdxj J ^ ^ dxidxj R ^ ^ 

The generalized derivative ([7]) can thus be seen as the mathematical operation using 
which the 2nd-order differentiation with respect to components of r of a function 
defined as the integral with respect to r' of an integrand involving the factor 
l/R = l/\r — r'\ may be performed under the integral sign: 
2 r „t*.i\ r a2 



9 2 f^P(r') [ A ,, f ,, d 2 1 



dxidxj J R J dxidxj R 1 

where the 'density' p(r) now plays the role of a well-behaved test function. 

Expression ([2]) is a special case r' = (so that R = i — r' = r) of the generalized 
derivative (0, 

B 2 1 3tt — r 2 S 47r 

wlim 1 3 ^ 13 8(r - e) - — 5{r)5 tr (10) 



dxidxj r 

+ As the informal proof in [3] , informal derivations of the full Laplacian relation JT} typically use 
the divergence theorem with a function that has a 1/r 2 singularity at r = 0; for a more rigorous 
alternative see |5]. 
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Adding expressions (|10|) with i = j = 1. 2, 3, the non-delta- function terms cancel out, 
and we obtain the full Laplacian ([TJ as 

\/ 2 - = -iirSfr), (11) 
r 

where 

dx\ dx\ dx\ 

is now the generalized Laplacian operator. 

The generalized derivative (|10j) is the weak limit a —¥ of the corresponding 
classical derivative of the potential 

M r) = Jd 3 r'^l, R = \r-r>\ (13) 

that is due to a localized density p a (r) the spatial extension of which depends on a 
parameter a so that 



lim / d 3 r/(r)p a (r)=/(0), (14) 

where /(r) is any well-behaved test function; this condition is transcribed formally as 
wlim p a {r) = S(r). (15) 

a— >0 

This can be shown by calculating the limit a —> of the integral of the product of © 
with <p(r) = <f>a(r) and a well-behaved function f(r). It suffices to consider the case 
i=j = 1: 

lim /d 3 r/(r)^^ 

= lim / d 3 r/(r) lim / dVp 8 (r') 3(xi ~^ )2 ~ fi2 e(iZ-e) - ^ lim / d 3 r/(r)p a (r) 
= limlim| d 3 r/(r) | dV p a (r') ^ ~ ^ = - 9(iZ - e) - y | d 3 r/(r)<5(r) 
= lim lim / dV Po (r') / d 3 r f(r) ^ - — Q(R - e) - ^ /(0) 

a-s-0 e->0 J J H° 3 

- lim / dV p a (r')F(r') - ^ /(0), (16) 

a— >0 _/ d 



F(r') = lim / d 3 rf(r) 3{X1 ^ ^ G(R-e). (17) 



where 

TPf-r'\ - lim f H 3 r- fC^^lZ 

The function F(r') is well-behaved, and thus 

limy" dVp„(r')F(r') = J dV S(r')F(r') 
= F(0) 

/o 2 2 
d 3 r/(r)^^ 9(r-e). (18) 

The algebraic manipulations in (|16|) are legitimate operations of moving limits without 
changing their sequential order and of interchanging the orders of integration. 
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Equations similar to (|16|) and (1181) obviously hold for all the derivatives 
d 2 cj) a (r)/dx i dxj, i.j — 1,2,3. We thus have that, for any well-behaved test function 

/(<■), 

which is expressed formally as 



,. d 2 4> a {r) .. ix l x 1 -r 2 8 ij 4tt 
m .'. -8x~dx- = W H — ^ 1 Hi '- " T ^ ' 

B 2 l 



dxidxj 



(20) 



Since the limit a —> of the potential a (r) itself is the potential 1/r of a point 
density S(r), this result is the formal underpinning of a natural interpretation of the 
generalized derivative (B 2 / dxidx / r) as the 2nd-order derivative of the potential of 
a point source. Note that while the regularization used in (|20j) is of the spherical kind, 
the extended density p a (r) that generates the potential (f) a (r) does not have to have 
any particular symmetry as a function of r. Using the powerful methods and results 
of the theory of generalized functions and derivatives, the result ([217)1 can be obtained 
almost immediately (see [12] , Appendix) , but the approach adopted here required only 
a little more effort. 



3. Non-spherical regularizations 

3.1. Spheroidal regularization 

Let us replace the spherical excluded integration domain used in ([5]), ([7]) and (fTU|) by 
a spheroid of semiaxes e/7, e and e along the X\, Xi and x% axes, respectively, with 
the parameter 7 given by 

J=^=L=, 0<H<1. (21) 
V 1 — ir 

Regularization using such a domain may be suitable in applications involving effects 
of special relativity, according to which a spherical charge of radius e contracts to 
an oblate spheroid of this geometry when it is set in motion with a speed vc along 
the x\ axis. To find the modification of expression (|10p for the generalized derivative 
(a 2 1 'dxidxj)(l J V) that such spheroidal regularization entails, we only need to evaluate 
the difference 

lim / d\f(r) ZX ' X3 ~ rHl3 1 - lim / d 3 r /(r) ^ T ^ 

= ^J u ^rf(r) 3X ^- r% \ (22) 

where f(r) is again a well-behaved test function and the integration domain U e is the 
region delimited by the oblate surface ^ 2 x\ + x\ + x\ = e 2 and the spherical surface 
r 2 = x\ + x\ + xj = e 2 : 

U t = {(x 1 ,x 2 ,x 3 )- 1 -f 2 x 2 1 +x\ + x\ > e 2 r\x 2 + xj+ xl < e 2 }. (23) 
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As e tends to zero, the integration domain U t gets progressively smaller and closer 
to the origin r = so that, for any r £ U e , f(r) —> /(0) as e —> 0. The right-hand 
side of (f2"2")) can therefore be written as 

lim / dV/ir) ^:^ = /(0)lim / d 3 r ^ T ^ . (24) 

Here, the integral on the right-hand side can be evaluated easily in spherical 
coordinates. With x\ axis as the polar axis and cos# = £, we obtain for i = j = 1: 

rb J-i y e /vi+(7 3 -i)c 2 r 



tt J d£(3£ 2 -l)ln[l + ( 7 2 -l)£ 2 ] 

„ / 2 2arcsinu 2\ . . 

27F 3 o • 25 



The case i = j = 2 gives 

„2 „2 i r 1 r 2 ^ 



j 3^_ I I ^ I ^ [3(1 _ ^ cog2 ^ _ x] ln[i + ( ^ _ m 
Ju e r 1 J-i Jo 



1 d£(l-3£ 2 )ln[l + ( 7 2 -lK 2 ] 



„ ( \ 1 arcsiniA . . 

= 27r U^ + ^J' (26) 

and the same result obviously will be obtained for i = j = 3. The mixed cases i =/= j 
will all yield zero on account of the integration with respect to the azimuthal angle 
<fi. A notable feature of the results ([25]) and (|26|) is that they are independent of e. 
Collecting all these results, equation (|24l) can be written as 

lim / d 3 r/(r) ^ = ^ = 2?^) - f%]/(0), (27) 

where 

!2/v 2 — (2/71; 3 ) arcsinw i = j = 1 

1 - 1/u 2 + (l/7v 3 )arcsini> i = j = 2,3 (28) 
i^j. 

This result establishes a generalized-function identity 
wlim 3XlX3 ' ~ ^ [6(7 2 x 2 + x 2 + x 2 - e 2 ) - 9(r - e)] = 27r^(«) - § %]*(r), (29) 

using which the difference (|22"T) with any well-behaved function f(r) can be evaluated 
immediately. According to this identity 

wlim "^ XlX i r ^ Qi r _ e \ 

e^O r 5 

= wlim 3 ^ -^ e(7 2 x 2 + * 2 + 4 e 2 ) - 27r[ ff<i («) - § %]<5(r) (30) 
and thus the generalized derivative (JTOJ) can be re-written as 

1 = wlim 3X ^ r 7 2 ^ 9(7 2 x 2 + x\ + x\ e 2 ) - 2ic gij (v)6(r), (31) 
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which is a spheroidal-regularization alternative to the standard, spherical- 
regularization expression (fTU|) for the generalized derivative (d 2 /dxidxj)(l/r); (f3"Tj) 
reduces to (|10[) in the limit v — > since lim^oT = 1 an d hm„_yo 9ij ( v ) — f <%■ An 
electrodynamic example in which expression pip is useful is given in Appendix. 



3.2. Cylindrical regularization 

In some applications, cylindrical coordinates are natural to the problem and the 
requisite integrations are performed most easily in these coordinates. We can find 
the expression for (8 /dxidxj)(l/r) that employs a cylindrical regularization by 
evaluating the difference 




d 3 r/(r) 



33? ? 3- 



1^3 



f d 3 r/(r 



3 3? ^ 3?^ T* & % j 



(32) 



Here, the integration domain Te, K is defined as 

Te, K = {(xi,x 2 ,x 3 ); \xi\ > neUxj + x\ > e 2 }, (33) 
which is the complement of a cylinder of base radius e and half-height ne, parallel to 
the X\ axis and centered at the origin r = 0, and the integration domains Vein an d 
Ve^K are defined as 

= {{x 1 ,x 2 ,x 3 );xj+x 2 +x 2 >e 2 nx 2 + x 2 <e 2 n\x 1 \< ne} (34) 
V e ( f2 = {[xi,x 2 ,x 3 );xl +xl + xl < e 2 C\ x 2 2 + xl < e 2 n\x x \ > kc}, (35) 
which are regions delimited by the surface of the cylinder and the spherical surface 
r X-^ ~\~ x^ x 3 — — e 2 . The region vif^ is nonempty only when n < 1, in which case 
the integral over vi 2 « has to be subtracted from that over the region Ve,« ■ 

As e tends to zero, the regions Ve,K progressively shrink and collapse onto the 

(1 2) 

origin r = so that, for any r £ Ve,« , /(r) — > /(0) as e — » 0, and thus 



lim 



d 3 r 



/(0) lim 



d 3 r 



fir 



3 3? ^ 3/ j 



d 3 r 



d 3 r 



33?^ X j 



(36) 



where the large brackets are used to denote the difference of the indicated integrals. 
The integrals on the right-hand side are evaluated easily in cylindrical coordinates 
s, cf), x\. The case i = j = 1 gives 



d 3 r 



d 3 r 



(2) 



3x 2 - 



4ir sds 
'o 



2x 2 - s 2 
sJT^^s 1 (x 2 + s 2 f/ 2 



= 47T 



The cases i = j = 2, 3 give 



d 3 r 



3.r 



2,3 



2 1 sds 
Jo JV^ 



dxi 



2tv 



(37) 

(3cos 2 0- l)s 2 -x\ 

( X 2 +s 2)5/2 
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i f ^ 23N 

= 2n / sas axi 



{x\ + s 2 ) 5 /2 



-*(§ (38 » 

and the mixed cases i ^ j yield zero because of the integration with respect to <fi. 

Similarly to the establishing of identity ([30|) . these results now establish a 
generalized-function identity 

3xiXj v 5{j . . 
whm — 0(r — e) 

= wlim 3xlX i 7 r ^ [0(bi| - ne) + Q{x\ + x 2 - e 2 )Q{ne - \ Xl \)\ 
e— >o r" 

-27r[ft ii («)-|5y]<y(r) ) (39) 
where 



2 - 2k/ vl + k 2 i=j = l 
hij(K) = { k/VI + k 2 i = j = 2, 3 (40) 

i ^ j. 

Using identity (|39[) . the generalized derivative (fT0|) can be written as 

/; " I = w i im ~ r2 ^" [6(1^1 - ne) + Q{x 2 + x 2 - e 2 )6(«e - \ Xl \)] 



dxidxj r e^o r 

-27rft ij -(/s) ( y(r). (41) 

This is a cylindrical-regularization alternative to (ITU1) . In the limit k — > 0, the delta- 
function term in (|4"T]) simplihes to — 47r5(r) for i = j = 1 ([J], p 29) and to otherwise; 
such regularization is implemented automatically by using the cylindrical coordinates 
s, 4>, x\ and performing the requisite integration over the whole space ]R 3 , but so that 
the integration with respect to the variable x\ is done last: 

dr f( r )-^—5 = / dx i[ sds d0/(s,0,xi) 



dxidxj r J_ 00 \J J ' {x\ + s 2 f/ 2 , 

-4^/(0)^1^1, (42) 

where dij (s, <fi, x\ ) is the function obtained by the transformation of 3x{Xj — r 2 6ij from 
the Cartesian to the cylindrical coordinates (e.g., du(s, (j), x%) = 2x\ — s 2 ). 



4. Concluding remarks 

In reference |12) . the equivalence of the spherical- and spheroidal-regularization 
expressions for the generalized 2nd-order partial derivatives of 1/r was illustrated 
by a relatively laborious explicit calculation of the weak limit a — ¥ of the derivatives 
of the Coulomb potential of a charged conducting spheroid of Unite extension a. The 
derivation of the generalized derivatives given here, together with the presented results 
on non-spherical regularization, ensures that any similar explicit calculation must yield 
the same result. 

A point worth making is that while the standard, spherical-regularization 
expression (fTOj) and the non-spherical-regularization expressions (f3Tj) and (|4"Tj) for the 
generalized derivatives (a 2 /dxidxj)(\/r) are guaranteed to yield the same results in 
an integral with a well-behaved function, care should be taken in numerical work to 
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use an integration grid that is compatible with the kind of regularization employed 

Very recently Frahm's formula @ has been criticized as being valid only when 
averaged over smooth functions, and, to remedy that, an expression 

d 2 1 = ZxiXj - r 2 5 l3 _ <5(r) (43) 



dxidxj 



has been proposed [BJ. However, our analysis shows that the only 'flaw' of Frahm's 
formula is that it does not indicate explicitly the spherical regularization that it 
assumes. We note that expression (|43)l still suffers from the lack of an appropriate 
regularization of the non-delta function term. Moreover, the delta-function term as it 
stands there is ill-defined; it would become meaningful in an integration in spherical 
coordinates and the replacement of the 3D delta function 8(r) by the radial equivalent 
8(r) / {Airr 2 ). Clearly, the 'general' expression (T4"3")) will yield correct results only in an 
integration in spherical coordinates, with the angular integration of the term involving 
the non-delta- function part done first. 
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Appendix 

We shall employ here a spheroidal-regularization generalized derivative in a calculation 
of the difference A^p — A^ between the xi-components of the Coulomb- and Lorenz- 
gauge vector potentials of a unit point charge moving uniformly with a velocity v 
along the x\ axis. 

The difference A&p — A^ satisfies an inhomogeneous wave equation [2J: 

□ [4? (r , t) - Ag> (r, t)} =~v^ (c) (r, t), (A.l) 

where □ = V 2 — d 2 /dt 2 is the d'Alembertian operator (we use Gaussian units with 
the speed of light c = 1) and 

0< c ) (r , t) = 1 (A.2) 

- vt) 2 +x%+ x% 

is the Coulomb-gauge scalar potential of the charge. But an inhomogeneous wave 
equation 

□/ = s{xi - vt,x 2 ,x 3 ), (A.3) 

whose source term is 'moving' with a constant velocity v along the x\ axis, can be 
simplified by a simple transformation of the variables to a Poisson equation, the 
vanishing-at-infinity solution of which is given in terms of the original variables by 

7 f A z, s{x' x -vt,x' 2 ,x'z) 



f(r,t) = -J- / dV "7 (A.4) 

4tt J Vrl^i - x i) + (as - 4) 2 + fa - x' 3 ) 2 
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where 7 = (1 — v 2 ) -1 / 2 [T3| (see also UUJ). The wave equation (|A.ip is therefore 
solved by 

47r J V7 (#1 - x\Y + (x 2 ~ x' 2 ) z + (x 3 - x' 3 j 2 

Note that, in order that the integral is defined properly, the 2nd-order derivative of 
the point-charge potential cb^ must be here a generalized one. However, instead of 
using the standard spherical-regularization expression for the derivative, it will be seen 
that the evaluation of the integral on the right-hand side of (|A.5|) is facilitated greatly 
when a sphcrodial-regularization expression is used. Following pip , this is given by 

"=w lim -f -*e<«--«) 

dx\ e->o R 5 y ' 

. ( 1 arcsiniA . . , . „. 

- 4tt — — S( Xl - vt)5(x 2 )S{x 3 ), A.6 

where 

R = J (xi - vt) 2 +xj+xl and R* = ^J-f 2 (xi - vt) 2 + x\ + x\. (A.7) 

The integration in ()A.5|) is now done in two steps. First, the term involving the 
delta- function term of (IA.6I) is integrated readily and a transformation 7(2^ —vt) —> x\ 
is performed in the remaining integrand. This gives 

4?(r,t) ~ AV(r,t) = G(r,t) - (l - ±, (A.8) 



where 



G(r, t) = lim / . - "V- ' ' ^ ^— ©(^ " e). (A.9) 

Here, ii* is a vector with components 7(2:1 — vt), £2, X3. The expansion of the factor 
1/l-R* — r'\ in Legendre polynomials now can be used to separate angular and radial 
integrations; moreover, the use of spherical coordinates r',6\(f>' with the angular 
integration done first implements the regularization limit e — > automatically. Thus 

V » ; 2i?* ^ 2/(2/ + 1) 7_ x ^ (1-^/2)5/2 M ^ h v 1 

where £* = 7(^1 — vt)/R*, £' = a^/V = cos#', and the result 
dr' r 21 4/ 4- 1 1 

— = + 1 _L z > 1, (A.11) 

r> r 2 ' +1 2/(2/ + 1) R* ~ y J 

where r< (r> ) is the lesser (greater) of r' and R* , is used. The summation in (|A.10[) 
runs only over Legendre polynomials of even non-zero order since the integration with 
respect to £' yields zero when the Legendre-polynomial order is zero or odd. 

Second, the independent part of the integrand is expanded in powers of v 2 to 
facilitate the integration with respect to 

( 3 - 2 ^ 2 - 1 = V ^# [(2, + 3)C' 2 - 2n - 1] C ' 2 V". (A.12) 
(l-w 2 £' 2 ) 5 /2 ^ (2n)\\ LV ;s js v 1 



2x' 1 2 / 7 2 - x 
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Using this expansion in (IA.10|) , we obtain after interchanging the orders of summation 
and integrating term by term with respect to 

«(M»™£[r(" + !)]V»'f < 4! + 1 »^"> 

V n=0 



^r(n-/ + 2)r(n + Z + |) 



1 oo 

— y 



,2 " 1)!!]2 [(2n + l)«*) 2 " ~v 2n ] 



vR* ^ (2n + l) 

n— 1 v ' 

arcsin v 




, „ . (A.13) 

y/l - v 2 ^ • 1 

Here, in the 1st line, the series over I terminates at I = n+1 since all its Z > n+l terms 
vanish; in the 2nd line, the terminated series is summed and the resulting series over n 
re-arranged so that it has an overall multiplier 1/v; and, in the 3rd line, the series over 
n is summed using the well-known expansions of the functions (l + x) -1 / 2 and arcsin a; 
in powers of x. Using (|A.13[) in (|A.8p and then the definition £* = 7(^1 — vt)/R* with 
the definition (|A.7[) of R*, we obtain finally 



A^(r,t)-A^(r,t) 



y/1 - V 2 £* 2 




y/(x! ~ Vtf + X 2 + X 2 ^J{ Xl -vt) 2 + (X 2 +X 2 3 )/y 



2 



(A.14) 



The same closed- form expression for the difference between the X\ components of 
the Coulomb- and Lorenz-gauge vector potentials of a point charge moving uniformly 
along the x\ axis was obtained in [2] by calculating the requisite gauge function for 
the transformation between the Lorenz and Coulomb gauges using a formula derived 
by Jackson [T5] . 
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